Valley-Hall Kink and Edge States in Multilayer Graphene 
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We report on a theoretical study of one-dimensional (ID) states localized at few-layer graphene system ribbon 
edges, and at interfaces between few-layer graphene systems with different valley Hall conductivities. These ID 
states are topologically protected when valley mixing is neglected. We address the influence on their properties 
of stacking arrangement, interface structure, and external electric field perpendicular to the layers. We find 
that ID states are generally absent at multilayer ribbon armchair direction edges, but present irrespective of 
crystallographic orientation at any internal valley-Hall interface of an ABC stacked multilayer. 



I. INTRODUCTION 

Metallic surface states in a system with an insulating bulk 
are often related to topological order. 1 An important example 
is provided by the quantum Hall effect of two-dimensional 
(2D) systems in presence of an external magnetic field, in 
which two-dimensional (ID) edge states 2 accompany integer 
valued Chern indices 3 of bulk 2D bands. The recent identifi- 
cation of topological insulator materials 4 has provided a new 
example, one which does not rely on external magnetic fields. 
In topological insulators strong spin-orbit interactions yield 
bulk bands that have non-trivial values of a Z2 topological in- 
dex and support topologically protected surface states. The 
present study explores a 2D bulk example in which the rel- 
evant topological index is, as in the quantum Hall case, an 
integer-valued Chern index even though no magnetic field is 
present. 

Our work is motivated by the suggestion of Martin, Blanter, 
and Morpurgc 5 that ID states can be induced in bilayer 
graphene by changing the sign of an inter-layer electric field. 
(Electric field below will always refert to a field directed be- 
tween the layers of a few-layer graphene system.) These states 
have a formal structure similar to that of the zero modes that 
appear in the A phase of 3 He thin films at domain walls be- 
tween regions with opposite spontaneous orbital moments. 6 
The ID states in Ref.[5] can be understood as being a con- 
sequence of separate Chern indices of opposite sign associ- 
ated with the K and K' Dirac points of bilayer graphene. The 
property that these Chern numbers are implied by the momen- 
tum space Berry curvature of bilayer graphene's bands 7 when 
the electric field is non-zero is referred to as the valley Hall 
effect. The valley Chern numbers emerge when the electric 
field breaks inversion symmetry to open a gap in the 2D bulk 
bilayer electronic structure. However, because the two val- 
leys share the Brillouin-zone of a bilayer graphene crystal, 
separate Chem numbers are never precisely defined and are 
strictly speaking an artifact of the commonly employed con- 
tinuum k ■ p electronic structure model. Correspondingly, the 
ID states that are the subject of Ref.[5] are not guaranteed to 
be present at all energies and are not topologically protected 
against perturbations that couple different valleys. The goal of 
the present study is to assess the degree to which these caveats 
are practically important. 

ID states of the type we consider were first found in numer- 
ical studies of zigzag bilayer ribbons 8 , and later recognized as 



valley Hall edge states by Morpurgo et al. 9 They have also 
been previously studied theoretically in single-layer graphene 
samples with an imposed staggered potential / ' 11 i.e. a po- 
tential that has opposite signs on the graphene honeycomb's 
A and B sublattices. A more detailed analysis of the cor- 
respondence between bulk and edge in bilayer graphene has 
been carried out recently, concluding that the existence of ID 
states at an edge depends on its morphology. 12 In the present 
paper we consider both bilayer graphene and other few layer 
graphene systems in which electric fields open up an energy 
gap and yield a valley Hall effect. Continuum model consid- 
erations suggest that ID states should be present at the edge of 
graphene ribbons with a valley Hall effect and at interfaces be- 
tween systems with different valley Hall effect quantum num- 
bers. We will refer to the former type of ID state as an edge 
state and to the latter as a kink state. In the case of bilayer 
graphene, for example, an interface that supports kink states 
is easily produced by changing the sign of the electric field 
along a line inside the material. To create the corresponding 
ID kink states in monolayer graphene, it would be necessary 
to change the sign of the staggered potential. 

In this paper we use 7T-band tight-binding models to assess 
the influence of stacking orders and edge geometries on edge 
and kink state properties. As in the bilayer case, we find that 
the valley Hall edge states in multilayers do not survive for 
armchair edge terminations. However, kink states are clearly 
present for both zigzag and armchair crystallographic orienta- 
tions of an internal interface along which the valley Hall quan- 
tum number changes. Section II contains the main results of 
our work. We start discussing the valley Hall conductivity in 
terms of the low energy continuum model of ABC stacked Al- 
layer graphene. We conclude that we can normally expect N 
ID kink state branches per valley localized along electric field 
sign-change lines. Then we use ft-orbital tight-binding cal- 
culations on multilayer ribbons to test the continuum model, 
presenting results for the energy bands of valley Hall edge and 
kink states for a number of different cases. Finally in Sec. Ill 
we close with a brief summary and a discussion of our find- 
ings. 



II. VALLEY-HALL EDGE AND KINK STATES IN 
MULTILAYER GRAPHENE 

The electronic structure results for multilayer graphene rib- 
bons presented here were obtained using a ^-orbital tight- 
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binding model Hamiltonian with nearest-neighbor hopping 
and a lattice position dependent external potential Uj 



a. 



H = L nJ c J c J ■ L u > 
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The hopping amplitude Yi.j we used is equal to t = 2.6 eV for 
in-plane hopping and t± — 0.34 eV for out of plane hopping. 
c] and cj are creation and destruction operators at i th and /'' 
lattice sites. This model can yield ID edge state or kink state 
branches at energies inside bulk gaps which can be explained 
qualitatively in terms of the model's bulk valley Hall effect. 
We start by discussing the valley Hall properties of multilayer 
systems using a low energy model Hamiltonian, before return- 
ing to tight-binding model calculations carried out for ABA 
and ABC stacked bilayers trilayers, and tetralayer graphene 
ribbons. We consider a variety of different examples that cap- 
ture some essential features of stacked multilayer electronic 
structure. 



A. Valley Hall effect and associated ID states in chirally 
stacked N-layer systems 

The low energy Hamiltonian in ABC stacked (rhombohe- 
dral) AMayer graphene is useful as a simplified route for gain- 
ing insight of the system despite the simplifying assumptions. 
In particular, the valley Hall properties of ABC stacked (rhom- 
bohedral) AMayer graphene can be derived from the low en- 
ergy band structure in which Bloch states are localized mainly 
on top and bottom layers. Because the phase difference be- 
tween top and bottom layer wavefunction components varies 
more rapidly with momentum direction (measured from the 
K and K' Dirac points) in larger N systems, the valley Chern 
numbers increase with Af.'^ilit To be more precise the wave 
functions of the states closest to the Fermi level reside mostly 
on the top and bottom layer lattice sites without a vertical 
neighbor, and the valley Chem number in the presence of an 
electric field is equal to N/2 except possibly at very low car- 
rier densities and weak electric fields where weak band struc- 
ture features can play a role. (Note that for odd N this quantity 
is not an integer; we nevertheless refer to the values as Chern 
indices for convenience.) The sites with vertical neighbors 
have more weight in higher energy bands. 

Strictly speaking, Chern numbers should be calculated by 
integrating Berry curvatures over the whole Brillouin zone, 
including contributions from near both K and K' points. This 
total Chern number always vanishes, as it must when time- 
reversal symmetry is not violated. In the continuum model 
approximation, however, we can speak of valley resolved con- 
tributions to Chern numbers because the valley indices are 
clearly distinguished. In this picture, familiar quantum Hall 
considerations suggest the presence of ID channels at edges 
and along lines where the valley Hall conductivity changes. 
The microscopic tight-binding model calculations we perform 
are intended to test the degree to which these considerations 
are reliable. 

Our working assumption in this section is that the following 
effective continuum model 13 captures essential properties of 
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FIG. 1: Schematic illustration of the relationship between valley 
Hall effects and ID conduction channels at interfaces expected on 
the basis of continuum model considerations. The number of ID 
modes per valley is an integer n evaluated from differences between 

valley Chern numbers n = A?| — . When the layer number A' is 

odd, the integral of the Berry curvatures is half-odd-integer, a prop- 
erty that is related to the half-quantized Hall effect of Dirac systems. 
We find that the number of ID channels is then usually reduced to the 
integer part of the integrated Berry curvature difference. The sense 
of the arrows directed perpendicular to the page indicates the sign of 
the Chern number associated with each valley. Upper panel: Illus- 
tration of ABC trilayer graphene with opposite electric field signs in 
left and right regions with three ID modes per valley. Lower panel: 
lunction formed between pentalayer and trilayer regions under a uni- 
form bias potential. In this case the valley Hall conductivities have 
the same sign but different magnitude on opposite sides of the inter- 
face and the number of ID channels is expected to be proportional to 
the difference of the individual Berry curvature integrals. The con- 
tinuum model picture illustrated here can be invalidated by atomic 
scale physics at the interface, particularly when the number of layers 
on opposite sides of the interface is different. 



the full Hamiltonian: 

Hn = —t± 
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(2) 



where — 1± is the interlayer hopping parameter, v = X>p%jt^, 
Vf ~ c/300 is the Fermi velocity of single layer graphene, 
% = h (k x + ikyty) where i z — ±1 labels the K and K' valleys, 
(k x ,k y ) is crystal momentum measured from a Dirac point, 
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and A is the potential difference between layers produced by 
the electric field. This Hamiltonian leads to momentum-space 
Berry curvatures that are sharply peaked near (k x ,k y ) = 0. 
The valley Chern numbers are obtained^^ by integrating 
the Berry curvature over 2D momenta (k x ,k y ) continued to 
°° to obtain C x . — A^T z sgn(A)/2. The valley Chern numbers 
are sometimes referred to as topological charges, and de- 
noted by N3 as a reminder that the topology indices of gapped 
2D systems may be thought of as a dimensional reduction 
of topological charge N3 at gapless 3D Fermi points^ The 
contribution to the Hall conductivity from a particular val- 
ley is (7§ = N$e 2 /h. In a continuum model, the number of 
one-dimensional (ID) channels per valley at an interface be- 
tween two bulk regions with different valley Hall conductiv- 
ities cr^ = N^e 2 /h is equal to the difference between their 
valley Chern numbers . 5 ' 6 ' 11 For ballistic transport each chan- 
nel contributes e 2 /h to the two-probe conductance. When the 
sign of the electric field is reversed at an interface, the Chern 
number difference is 2 x (N/2 )= N, equal to the layer number 
as illustrated schematically in Fig. 1. 

The generalization of the notion of a valley Hall conduc- 
tivities from one and two layer systems to general N layer 
systems can also be made using the explicit equations for the 
ID interface states by following a procedure similar to that 
outlined in Ref . j5[ : 

-V(x)u + K N (d x + k y ) N v = eu 
K N (d x - kyf u + V(x)v = ev. (3) 

Here — —t± (~iv>Fh/t±) N , N is the number of layers in 
the system, and V(x) is a general position-dependent function 
which specifies the difference between top and bottom layer 
potentials. In the following subsections, we will present tight- 
binding calculations for multilayers for a variety of different 
external potential profiles and discuss the validity of the qual- 
itative picture summarized in Fig. 1. 

B. Multilayer ribbons under a uniform electric field 

The simplest example of valley Hall edge states are those 
that appear in ABC stacked multilayers under a uniform ex- 
ternal electric field. Tight-binding band structure of bilayer 
graphene ribbons by Castro et al.&lS- has demonstrated the 
presence of metallic edge states which cross the Fermi level 
in neutral zigzag terminated bilayer ribbons and in chirally 
stacked multilayer zigzag ribbons. The number of valley Hall 
edge states branches in each propagation direction is equal 
to [A/2] 15 . Changing either edge termination or stacking 
sequence introduces qualitative changes in the ribbon band 
structure as we now describe. 

In Fig. 2, we plot band structures for ABC and ABA zigzag 
and armchair terminated trilayer and tetralayer ribbons. In 
all these calculations we have maintained the same interlayer 
potential difference A/t = 0.1 between top and bottom lay- 
ers and have taken the electric field to be uniform. In the 
tetralayer case the smaller potential difference A' — A/3 be- 
tween the inner layers to facilitate distinguishing their edge 
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FIG. 2: Band structures of ABC and ABA stacked zigzag and arm- 
chair terminated trilayer and tetralayer ribbons in the presence of a 
uniform electric field, (a) ABC stacked layers with zigzag edges; 
(b) ABA stacked layers with zigzag edges; (c) ABC stacked lay- 
ers with armchair edges and (d) ABA stacked layers with armchair 
edges. In each subfigure, the left and right panels represent trilayer 
and tetralayer ribbons, respectively. For ABC stacking the electric 
field opens a bulk band gap containing edge states in the zigzag case 
but not in the armchair case. The number of edge state branches is 
independent of the ribbon width, whereas the number of bulk state 
branches is proportional to ribbon width. In the case of ABA stack- 
ing, ribbons are metallic in the trilayer whereas a small gap opens in 
the tetralayer geometry. 



state bands associated with respect to the outer ones. In mul- 
tilayer zigzag ribbons with ABC stacking, the uniform elec- 
tric field generates a band gap in the sample bulk 14 with edge 
states in the gaps. In agreement with previous analysis^ 5 -, we 
find [N/2] valley Hall edge states in an A-layer ribbon in each 
valley that propagate in opposite directions and are localized 
at opposite edges. Hence, for zigzag edges the number of ID 
states at the physical boundaries of the ribbons with vacuum 
can also be discussed in terms of valley Chem number dif- 
ferences between the bulk region and the vacuum N™ c = 0. 
In the armchair termination both K and K' valleys appear at 
the same projected ID momentum and the Hall edge states 
are annihilated. This observation is helpful in distinguishing 
edge and kink states in ribbons with internal electric field sign 
changes. 

For biased ABA multilayer zigzag ribbons the system also 
has conducting edge states, but their pattern is more complex 
than in the ABC case. (Bulk and edge ribbon bands can be 
distinguished by their dependence on ribbon width.) In ABA 
stacked trilayer ribbons an external bias increases the num- 
ber of bulk channels, while it opens a bulk gap in the ABC 
case^il 
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C. ABC stacked zigzag and armchair multilayer ribbons with 
kink states 



We have just seen that a biased bilayer and ABC trilayer 
graphene ribbons are gapped in the bulk but have metallic val- 
ley Hall edge states crossing the Fermi level when they have 
a zigzag edge termination. Now we consider ribbons with a 
electric field sign change at the ribbon center that is expected 
to produce AMD channels in each valley. In these ribbons, 
the regions near the ribbon border will show valley Hall edge 
states similar to those present in a bilayer with a uniform elec- 
tric field. In Fig. 3 we plot the band structures of bilayer, 
trilayer and tetralayer ribbons with zigzag edge terminations 
subject to a step-like interlayer potential discussed earlier. 

In zigzag terminated ribbon geometries the K and K' mo- 
mentum projections appear at the two valley points located at 
k = 2n/3a,4n /3a and can therefore easily be distinguished. 18 
In agreement with the continuum model analysis presented in 
Sec. II, we see N bands of confined kink states with a well 
defined propagation direction for each valley K or K' crossing 
the Fermi level. Each valley has doubly degenerate additional 
metallic edge state branches, with a velocity opposite to that of 
the confined states. In the uniformly biased case, edge states 
in a given valley that are localized on opposite edges have op- 
posite propagation directions, whereas they propagate in the 
same direction when a kink is present. For ribbons with inver- 
sion symmetry at the ribbon center, the co-propagating edge 
state channels are degenerate. 

As mentioned previously, the projection of the 2D bands 
of graphene to obtain the ribbon band structure places K and 
K' valleys at the same momentum in the armchair edge case. 
Therefore, unlike the case of the zigzag ribbons, it is not pos- 
sible to identify valley labels from ribbon band structure plots. 
In the case of edge states, these difference eliminates the ID 
channels completely. As we see in in Fig. 3 this is not the 
case for ribbon states which appear to be as robust in arm- 
chair and zigzag directions. When the electric field profile has 
sharp spatial variation, there is a barely visible gap opening, 
similar to the one found in direction of graphene under a stag- 
gering potential 10 . This gap size decreases quickly when the 
potential variation at the domain becomes smoother. 

In the band structures of bilayer, trilayer and tetralayer 
graphene armchair ribbons subject to a kink step bias around 
the ribbon center, one can clearly identify two, three and four 
ID states for each propagation direction, corresponding to 
confined states at the domain wall for wave vectors around 
k ~ where both right and left going states coexist with simi- 
lar Bloch function wave vectors. In armchair edges we do not 
find edge localized states crossing the Fermi level as we had 
found for the zigzag terminated systems. 



m. SUMMARY AND CONCLUSIONS 

External electric fields between layers give rise to gaps 
at the carrier neutrality point in bilayer and ABC multilayer 
graphene systems. We have shown that ID transport channels 
appear along lines where the sign of the inversion symmetry 
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FIG. 3: (Color online) Upper Panel: Band structure of bilayer, tri- 
layer and tetralayer graphene zigzag ribbons with an electric field 
sign change at the ribbon center. We can clearly observe two, three 
and four ID valley Hall kink states with a common sign of velocity 
located around the valley point. The branches which correspond to 
wavefunctions localized at the ribbon center are plotted in red. The 
other branches that cross the bulk gap are doubly-degenerate edge 
states. Lower Panel: Band structure of bilayer, trilayer and tetralayer 
graphene armchair ribbons with an electric field sign change at the 
ribbon center. As in the zigzag case we can clearly identify two, three 
and four three ID kink state branches. In armchair edges we do not 
find ID edge state channels. 



breaking potential changes. This finding generalizes results 
obtained previously for the bilayer 5 and monolayer cases i 10 ! 11 
The number of these metallic ID kink state branches is pro- 
portional to the number of layers N and can be related to the 
bulk valley Hall conductivity. 15 . The states we have consid- 
ered arise at the boundary between two regions with oppo- 
site valley Hall conductivity. Because the Hall conductivity 
changes sign in opposite directions in the two valleys both 
produce ID kink states and they propagate in opposite direc- 
tions. 

Similar valley Hall effect considerations suggest that kink 
states should occur at boundaries between ribbons with dif- 
ferent thicknesses. Our % -orbital tight-binding model cal- 
culations for ribbons with a trilayer/tetralayer boundary find 
that the interface electronic structure depends on whether the 
stacking sequences is ABC or ABA. Bulk valley Hall effect 
values are also unreliable at the edge; in particular we find that 
edge states are absent in multilayer armchair terminated rib- 
bons, as found earlier in the bilayer or monolayer case. This 
finding is perhaps expected since the ID momentum projec- 
tion does not distinguish valleys in this case. It is the robust- 
ness of the kink states at internal electric field sign changes 
that is perhaps the surprise. It is important to determine if it 
persists in the presence of disorder and turns in the propaga- 
tion path. 
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